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The six elastic constants (and six clastic compliances) of corundum were determined in 
the kilocycle per second frequency range by an accurate resonance method. The results 
were checked in the megacycle per second range with a less accurate, pulse velocity method. 
The elastic moduli for polycrystalline alumina calculated from the single crystal compliances 
determined by the resonance method are in good agreement with experimental values ob- 
tained on high density polycrystalline alumina. The variation of Young's modulus and of 
the shear modulus with orientation was calculated from the compliances and the results 
are shown graphically. The results of the present work do not agree well with previous 
work on single crystal sapphire. The specification of orientation and the theory used to 
calculate the elastic constants are given in detail to support the contention that the results 
of the present work are correct. 



1. Introduction 

The purpose of the present paper is to present a 
set of values for the elastic constants of single crystal 
corundum 1 at room temperature and to argue that 
these values are correct despite the fact that they 
disagree with values reported in three other investi- 
gations. Fifty-six independent measurements on 
29 different specimens yielded consistent results for 
the six elastic constants which this material should 
have. 2 

The complete set of six elastic compliances and 
six elastic constants were first reported by Sum lain 
Rao [l] 3 for synthetic sapphire. A similar deter- 
mination was made by Bhimasenachar [2] using 
natural sapphire; very good agreement with Sundara 
Rao was found except for c 33 . More recently, Mayer 
and Hiedemann [3] have redetermined the elastic 
constants and compliances for synthetic sapphire 
with results which are rather different from the 
previous investigators' values. Mayer and Hiede- 
mann suggest possible causes of incorrect resonance 
frequency measurements in the experimental method 
used by the Indian workers. 

The present work began with measurements of 
Young's modulus which were intended as the basis 
of a study of the changes in elastic moduli and 
internal friction as a function of temperature. These 
Young's modulus values permit calculation of four 
compliance constants which were found to be incon- 
sistent with either of the previously reported sets of 
compliance constants. All of the previous work was 



i Corundum is pure aluminum oxide, AI2O3, in the stable, a phase. The term 
sapphire is used synonymously in this paper. 

- A generalized theory, due to Laval and Raman, requires 10 constants and is 
discussed later. 

3 Figures in brackets indicate the literature references at the end of this paper. 



done at frequencies in the megacycle per second range, 
but the present determination of Young's modulus 
values was made in the low kilocycle per second 
range. A possible frequency dependence thus re- 
quired consideration. The present authors under- 
took a broadened program with the objectives of 
determining a complete set of elastic constants in the 
kilocycle per second range and checking the results 
in the megacycle per second range. The results of 
the broadened investigation are reported in this 
paper. 

In all the work discussed so far the classical theory 
of elasticity was used. This classical theory has 
been considered correct for more than 50 years and 
is still accepted by most scientists. Briefly, it 
assumes a linear relationship (Hookc's law) between 
6 independent components of stress and 6 independ- 
ent components of strain. The theory leads to 21 
independent elastic constants for the most general 
nonisotropic medium such as a crystal with no 
symmetry, and to 6 independent constants for a 
crystal with the symmetry of sapphire. This theory 
was challenged in 1951 by Laval who asserted that 
there are 9 independent components of stress and 9 
independent components of strain which would lead 
to 45 independent elastic constants for a crystal 
with no symmetry. This idea was developed by 
Laval and Le Corre, and independently by Raman 
and his collaborators. A good list of references is 
given by Joel and Wooster [4]. For the case of 
crystals with the symmetry of sapphire, Raman [5] 
asserts that there are 10 independent constants 
instead of 6 as in the classical theory. If this gen- 
eralized theory were correct, discrepancies in reported 
values of elastic constants might result from analysis 
of data in terms of a 6 constant theory when a 10- 
constant theory should have been used. However, 
results of the present work show that the classical, 
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6-constant theory represents the elastic behavior of 
sapphire in the kilocycle frequency range quite well. 
All equations in this paper are based on the classical 
theory. Because the present results do disagree 
with previous work, this paper presents the theory 
and experimental procedure in some detail. 

2. Specification of Crystal Orientation 

The values of the elastic constants (or compliances) 
of a nonisotropic material depend upon the choice 
of coordinate system. The use of a rectangular 
coordinate system is conventional. For each crystal- 
line material, the orientation of this rectangular 
coordinate system is usually chosen to take maximum 
advantage of its point group symmetry, and the use 
of this system is implied when elastic constants are 
discussed. In this section the coordinate system 
appropriate for sapphire (called the x^x-yXs system) is 
described in terms of the point group symmetry and 
the back reflection Laue pattern. Elastic property 
measurements on a crystal must be related to this 
coordinate system in order to interpret the results 
of the measurements in terms of elastic constants. 
The line in a specimen along which Young's modulus, 
the shear modulus, or the velocity of sound was 
measured, will be referred to as the specimen axis 
(SA). It is convenient to distinguish between the 
opposite ends of this line and the reference end will 
be taken as positive ( + SA). 

The point group of sapphire is 3 m [6]. Figure 1 
is a stereographic projection showing the symmetry 
elements, the rectangular system XiX 2 x s to which the 
elastic constants are referred, and the hexagonal 
coordinate system xyuz which is used for Miller- 
Bra vais indices. This notation for both systems is 
used by Nye [7] and will be used throughout this 
paper to refer to these coordinate systems. 

Figure 2 shows a stereographic projection of 
prominent poles in a back reflection Laue pattern of 
sapphire. Figure 2a shows the letter symbols of 




Figure 1. Stereogram of point group 3 m. 

The heavy lines indicate mirror planes, the open triangle indicates the '6 axis, 
and the diad symbols indicate 2-fold axes. 3:12:22:3 is the right-handed coordinate 
system used for elastic constant specification, x y u z is the Miller-Bra vais 
coordinate system. 

morphological crystallography [8] which can be 
assigned by measuring the angles between poles. 
This process is without ambiguity and is described in 
section 6. Figure 2b shows the Miller-Bravais 
indices for the conventional choice of coordinate 
system. Figure 2a is based on angular values ob- 
tained from Winchell [8], The recent X-ray pattern 
by Swanson and Fuyat [9] gives e/a= 1.365 for the 
morphological unit cell. The distinction between 
the morphological and the structural cells is discussed 
by Kronberg [10]. 

The XiX 2 x^ system is not uniquely defined by figure 
1 alone. Any 1 of 6 directions might be chosen for 
-\-Xii i.e., -\-Xi might be chosen in either direction 
along any of the three 2-fold axes. However, 
opposite ends of a 2-fold axis can be distinguished 
even though point group 3m has a center of sym- 
metry. The opposite ends are designated -\-a and 

— a in figure 2a. The distinction between -\-a and 

— a is illustrated in figure 3. The three -\-a direc- 
tions can be identified by examining the neighboring 
points on the stereographic projection. The choice 
of which is to be called -\-a and which —a is, of 




Figure 2. Stereograms of poles of prominent planes in sapphire. 

Figure 2a shows the letter symbols which can be assigned by measuring the angles between poles. Figure 2b gives the Miller-Bravais indices which require i 
choice of coordinate system. 
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course, a convention. This choice is important 
because the sign of the elastic compliance <v 14 de- 
pends upon it. 4 No convention for this choice for a 
nonpiezoelectric crystal is ejven bv Nye [7] or by 
the Standards on Piezoelectric Crystals [11]. The 
convention used here has been chosen to agree with 
that implied by figure 151 from Phillips [12], by 
figure 2 from Kjonberg [10], and by figure 5 and 
table 1 from WInchell [8]. The angles specifying 
the orientation of the specimen axis are shown in 
figure 4. Theta, the eolatitudo, or zenith angle, is 
the angle between -4 -x z and SA. Phi is the angle 



in the ./,./•. plane from + ./', to the projected specimen 
axis (PSA). The angle 6 is exactly the same as the 
angle p used by Winehell [8]. However, his angle <j> 
is measured from —a instead of from -f-a and is 
measured clockwise (looking in the — x z direction). 
WinchelFs definition of is not used here because 
we wish to adhere to the usual spherical polar angles 
in a right-handed system. 



'This can be seen by examining a coordinate system yiy^fz with y\ = —x\. If 
2/3=2-3, then one must have //?=— x% to maintain right-handedness. Examination 
of the transformation of aw shows that it is equal in magnitude but opposite in 
sign when expressed in the VxVtiji system. 




Figure 3. Distinction between +a and — a. 

The stereogram In figure 3a is centered on +a and shows thai the sequence of poles m trgc R 8 in the nearer minor planes has counterclockwise sense. Figure 
3b, centered on — a, shows t bat I he corresponding sequence has clockwise sense. 



1*3 




-a 



-f-a 




-+ a 



Figure 4. The angles used in specifying the direct ion of the specimen axis, SA, and its projection, PSA , in the XiXa plane. 

4a shows a perspective Of the important orientation angles. 4b shOWSa stcreographic representation of the orientation angles. 
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In discussing the range of d and <j> it is convenient 
to introduce the concept of the asymmetric triangle, 
which is the largest spherical triangle on a stereogram 
not containing any symmetry related directions. 
Any direction outside this triangle has a symmetry 
related direction within it. 

The range O°<0<9O° and -30° <0<3O° is suffi- 
cient to specify the orientation. That is, on a 
stereographic projection the range O°<0<9O° and 

— 3O°<0<3O° defines an asymmetric spherical tri- 
angle. However, this range of <f> is not sufficient 
to distinguish between equivalent points with mirror- 
Image environment. For example, the environment 
of (1120) is the mirror image of the environment of 
(2110). To distinguish such directions a range of 
O°<0<9O° and -30°<<£<90° is useful, and Kebler 
and Rudness [13] have recommended the use of this 
range in specifying the orientation. The range of 

— 30° <</>< +30° is adequate for specifying any 
physical property, such as elastic properties, which 
depend only on the direction of measurement and 
not upon the left-handed or right-handed distribu- 
tion of neighboring points. 

In carrying out the actual measurement of </> either 
of two equivalent methods may be used. In the 
first method <j> is measured from PSA to x x and the 
fact that 0=30° + 8 is equivalent to <£=30° — 8 is 
used to reduce any values larger than 30° to the 
range — 30°<<£<30°. In the second method # is 
measured from PSA to the nearest a-type direction 
without regard to whether it is +a or —a. The 
rule shown in figure 5 is then^used to determine the 
sign of </>. 




PSA 



- 0s o^ • 




Figure 5. Stereograms illustrating the sign convention for <f>. 

Note that the position of SA must be considered. The sign of <£ cannot be 
determined from the position of PSA alone. 



3. Relation of Young's Modulus and Shear 
Modulus to Elastic Compliances and to 
Orientation 

The stresses a t , and the strains, e iy will be written 
in Nye's matrix notation [7]. Hooke's law can be 
written in terms of the elastic constants, c ih as 



0~f ZrCij6j 

J 



(1) 



or in terms of the elastic compliances, s#, as 

e i =Xs ij (r j (2) 

j 

where the matrix of the % is the inverse of the 
matrix of the c^. For any material, c^=c^ (and 
Sij=Sji) which reduces the number of independent 
constants (or compliances) to 21. For sapphire, 
the symmetry of point group 3m further reduces 
the number of independent elastic constants (or 
compliances) to six and the resulting matrices are 
shown in table 1. 

Table 1. Elasticity matrices 



General elastic constant matrix 


en 


CIS 


Cl3 


Cl4 Cl5 


C16 


Cl2 


C22 


C23 


C24 C25 


C26 


Ci3 


C23 


C33 


C34 C35 


cm 


Cu 


C24 


C34 


C44 C45 


Cm 


cu 


C25 


C35 


C45 C55 


C56 


Cu 


C26 


C36 


C46 C56 


C66 


Elastic constant matrix for sapphire referred to the 






coordinate system 1112X3 




Cu 


C\% 


Cl3 


CH 





C\% 


C\\ 


C13 


-C14 





Cl3 


cu 


C33 








CU 


—Cu 





C44 














C44 


cu 











cu 


M (C11-C12) 


Elastic 


com pi 


iance matrix for sapphire 


referred to 




the coordinate system 2:1122:3 




5il 


su 


513 


su 





5l2 


811 


8a 


-su 





Sl3 


$13 


S33 








5l4 


—su 





S44 














$44 


2«14 











2514 


2(5ii— S12) 



In relating the elastic compliances to Young's 
modulus, it is necessary to distinguish between 
the "free" Young's modulus, E f , and the "pure" 
Young's modulus, E p , when analyzing flexural or 
torsional tests. The free Young's modulus is the 
value obtained when the specimen is completely 
free to deform elastically under the applied tensile 
stress. The pure Young's modulus is the value 
obtained when the specimen is tested in flexure and 
is prevented from twisting. In an isotropic medium 
the free and pure moduli are identical. Calculation 
of Young's modulus from flexural vibrations of long, 
thin rods of nonisotropic material corresponds to 
measurement of E f . A similar distinction must 
be made with the shear modulus, G. For rods of 
nonisotropic material of the dimensions used in this 
work, the modulus determined from torsional vibra- 
tions is very accurately the pure shear modulus, G v . 
The proof that E f and G p are the measured quantities 
involves the substitution in equations given by 
Hearmon [14] and by Brown [15]. We shall omit 
the proof here and give the results connecting the 
quantities determined experimentally, E f and G V) 
with the elastic compliances. Consider a rectangular 
coordinate system x[x'^x' z with the x' s axis along 
SA. For stress application along the x' z axis 
[14, 15]: 
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/■■■ 



s \\ I s :,:, 



2 



where 



"g; 

..'21 ..'2 
_ «Ss(«44+««)' 



(3) 

(4) 
(5) 

(6) 



Equations (4), (5), arid (6) can be combined to give 



'2_lo'2 



G f G v 2s 33 



which is a more useful relation between G lt and 
than eq (5). 



(7) 
G, 



The primed quantities in eqs (3), (4), and (7) 
must be related to the orientation and to the elastic 
compliances measured in the x x x 2 x z system. It can 
be shown [14] that the values of s m Su+s'm, and 
s'zI+s'm depend only upon the direction of z£ and not 
upon the direction of x[ and x 2 . Accordingly x\ and 
X 2 were chosen to simplify the transformation equa- 
tions as much as possible subject to the condition 
thai x[x 2 Xz form a right-handed rectangular set. This 
was accomplished by choosing x 2 to lie in the x x x 2 
plane, then choosing ,r\ orthogonal to x' 2 and x' z . The 
relation of x[x 2 x f z to x x x 2 x z is shown in figure 6. 

Equation 25 on page 70 and equation 10 on page 61 
of Cady [16] give the general transformation for s£ 3 
and #44 + .V55 respectively. An error in sign was found 
in 12 of the last 18 terms of the equation for #44+^5. 

esulting expressions 



When this error is corrected, the 
a re : 

1 



-E,(0,4>) 



s u sin 4 0+s 33 cos 4 



+ (2s 13 +s 44 ) sin 2 0cos 2 



+ 2.s 14 sin 3 cos0 sin 3$, 



(8) 



and 



^44 + ^55 



1 , , . 2flJ /1+COS 2 0\ 

+2 («u+«»-2«u-«44) sin 2 0cos 2 
— 4,s 14 sin 3 eos0 sin3</>, 



(9) 



which agree with eqs 48 and 49, page 76 of ( )ady [16] 
for trigonal symmetry. 

The transformations for s 34 and s 35 were obtained 
by using the general transformation law for fourth 
rank tensors, then reducing this to the trigonal 
symmetry of sapphire and making use of the particu- 
lar choice of x[x 2 described above. The resulting 



expressions are 

$34=3$i4 sin 2 cos 6 cos3<£, (10) 

and 

S35=2(2s 18 +$44— s u — 5 88 ) sin0 COS 3 

+ 2\s n -s n -^\ sin0 cos0 

+s u sin0 sin30 sin3<£. (11) 

A set of elastic compliances can be determined from 
these equations in the following manner: It is neces- 
sary to determine the values of E f and G p on at least 
four rods of known orientation. Equation (8) is 
written for each rod and the resulting system of 
simultaneous linear equations is solved for s n , s 33 , 
2.S13 + S44, and ,s I4 . These values are substituted into 
eqs (8), (10), and (11) and the results are used in 
eq (7) to calculate G f for each rod. These values of 
G f , 0, and 4> are substituted into eq (9) for each rod. 
The resulting system of simultaneous linear equa- 
tions is solved for 8n— s 12 , s 44 , Sn+S 33 — 2s 13 — s 44 , and 
s H . 

4. Relation of Elastic Moduli to Resonance 
Frequencies 

The values of E s and G p needed for the calculation 
of the elastic compliances were obtained by resonance 
frequency measurements on thin, circular rods. 
Young's modulus was calculated from the longi- 
tudinal resonance frequency and independently from 
the (lexural resonance frequency. The shear modulus 
was calculated from the torsional resonance fre- 
quency. For both of the Young's modulus calcula- 
tions, the equation relating the resonance frequency 
to the appropriate modulus is approximate, but the 
approximations are very good. 

For longitudinal vibrations, Young's modulus of 
an isotropic medium can be calculated from Ray- 
leigh's equation [17] which can be written 



(•> 2y,2\ 2 



(12) 



where pis the density in grams per cubic centimeter, 
/ is the length in centimeters, / is the longitudinal 
resonance frequency in cycles per second, a is 
Poisson's ratio, r is the radius in centimeters, and E 
is Young's modulus in dynes per square centimeter; 
1 dyne/cm 2 =10- 9 kilobars=1.020XlCT 6 kg/cm 2 = 
1.450X10 -5 lb/in. 2 . The term in parenthesis is the 
Rayleigh correction term for the finite thickness of 
the rod and neglects higher powers of r/l. Bancroft 
[18] has made an accurate calculation of the velocity 
of longitudinal waves as a function of r/X, where X is 
the wavelength which is equal to 2/ for resonance in 
the fundamental mode. The smallest value for 
which Bancroft gives a numerical result corresponds 
to r// = 0.05. Assuming Poisson's ratio is 0.25, the 
Rayleigh correction for this value is 1.00077 and 
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PSA m 





X' 

* 1 


X' 

A 2 


X* 
* 3 


*l 


cos -e- cos 


- sin $ 


sin cos 


x 2 


cos ■©■ sin 


cos $ 


sin •© sin 


X 3 


-sin -6- 





cos-6- 



Figure 6a is a 
direction cosines. 



Figure 6. The rotated coordinate system X1X2X3 and the direction cosine scheme. 

Figure 6b shows the relation stereographically. Figure 6c gives the values of the 



Bancroft's results give 1.00078. Evidently the 
Rayleigh correction is very good for small values of 
rjl. The largest value of rjl used in the present work 
was 0.0125 so that eq (12) is thought to be very 
accurate. For the range of rjl used in this work, the 
entire correction factor can be taken to be 1.0000. 
This equation was derived on the assumption of 



elastic isotropy, but the actual specimens were an- 
isotropic. Equation (12) gives E f for crystalline 
rods of the dimensions used. The value of Poisson's 
ratio varies somewhat for single crystals but for 
polycrystalline alumina is about 0.25, and since 
this only enters in the correction term, it seems that 
eq (12) can be used with this value of Poisson's 
ratio. 
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For flexural vibrations of an isotropic medium 
the best existing theory seems to be that based on 
an approximate 4 differential equation derived by 
Timoshenko [19] an 
Pickett [21]. Pickc 

of a correction Factor, 7', which multiplies tnc 
For an infinitely thin rod. His results can 
pressed as 



id studied by (locus 
's results are given 
tor, 1\ which multiplies 



[20] and by 

in t he Form 



result 
>e ex- 



#=0.31547 ^r, 



(13) 



where E is in dynes per square centimeter, and the 
symbols have the same meaning as in eq (12) except 
J which here is the flexural resonance frequency. 
In his original paper, Pickett [21] gives equations for 
calculating T For Poisson's ratio values of zero, 
1/6, and 1/3. A subsequent ASTM Standard [22] 
gives an interpolation Formula without giving the 
derivation or stating the range of validity. It 
seems best to use Pickett's original formula and 
interpolate using a quadratic approximation. The 
question is not significant for rods oF the size used in 
this investigation, the thinnest rods having r = 0. 050 
in. (r/Z=0.0083) and T= 1.0015, the thickest rods 
having r = 0.075 in. (r/Z=0.0125) and r=1.0033. 

The equation relating the shear modulus to the 
resonance frequency of an isotropic cylinder can be 
derived by considering the propagation oF torsional 
waves and requiring that the wavelength have a 
value such that standing waves are formed. The 
result is 



G=4 P IT, 



(14) 



where G is 
cent imeter, 
p is the density 
/ is the length. 



the shear modulus in dynes per square 
f is the torsional resonance frequency, 

n grams per square centimeter, and 
The significance of this equation for 



crystalline cylinders lias been considered by Brown 
[15] and by Ilearmon [14]. The details are compli- 
cated, hut the result is, as previously stated, thai for 
thin rods eq (14) gives Q v and eq (13) gives E f . 



5. Relation of Elastic Constants to the Veloc- 
ity of Ultrasonic Waves 

The theory required for the determination of the 
elastic constants of sapphire from resonance fre- 
quency measurement was given in sections 3 and 4. 
The present section presents the theory required for 
elastic constant determinations from wave velocity 
measurements. The velocity of sound as a function 
of specimen size has been studied by Tu, Brennan, 
and Sauer [23]. They Found that in cylindrical rods 
the velocity of longitudinal waves depended upon 
r/X where r is tin 4 radius of the rod and X is the wave- 
length of the sound wave. For values of r/X larger 
than 2.5, the longitudinal waves traveled at the speed 
of plane waves in an infinite medium. 

Similar results should hold for waves propagated 
in a rectangular block. The longest wavelength 
used in the present work was about 0.11-cm and the 
smallest block was 1.25-cm thick. These values give 



a thickness to wavelength ratio of 11.4. It seems 
that the theory for plane waves in an infinite medium 
can therefore be used. 

This theory is summarized by Kolsky [24] and by 
Mark ham [25]. The results are contained in a cubic 
equation in the variable x = pv 2 where p is the density 
and v is the velocity. This equation can conven- 
iently he written in determinant form 



(15) 



— X 


11 


G 


11 li-x 


F 


G 


¥ 


C-x 



= cnJ 2 -} - room- I C55W 2 \-2cidm-{'2c M mn+2ci5nl ) 

= <W 2 + C22M 2 + C44W 2 + 2c>^m -\- 2ci\inn -f- 2c^nl, 

z C55I 2 + cum 2 -f enn 2 -f 2crJm -\r 2c 3 4 run + 2c&n I , 

= c 56 / 2 + cum 2 -f c- M n 2 + (c 2 s -f c^) Im -f (c 23 + c u ) m n 

-\-(cz6+Ci&)nl, 

=Cul 2 +C4m 2 +C3 6 n 2 + {cu + cw)lm-\- (c 3 e+ C45) m n 

+ (Ci3+CM)nl, 



and 
H= 



C\d 2 ^C2&m 2 -^CA 5 n 2 -\r(c V 2 + CQ G ){m J r(c 2 5-{-Cu) ran 

+ (cu+Cn)nl. 



(16) 
(17) 
(18) 

(10) 

(20) 

(21) 



In these equations the c tj are the elastic constants 
and /, m, n arc 4 the direction cosines of the normal 
to the wave front. The distinction between velocity 
surface and wave surface for elastic waves has been 
discussed by Musgrave [26] and Markham [25]. It 
is sufficient to state that the measured velocities in 
the direction defined by Imn are those given by eq 
(15). It is important to notice that both the elastic 
constants and the direction cosines are referred to 
the same coordinate system. For sapphire, if the 
coordinate system is chosen to be the x x x 2 Xz system 
previously described, the matrix of elastic constants 
given in table 1 may be substituted to simplify eqs 
(16) through (21). This simplification has been done 
by Bhimasenachar [27] for calcite which has the same 
point group and consequently the same elastic con- 
stant matrix as sapphire. The resulting equations 
can be expressed directly in terms of the orientation 
by substituting £=sin0cos#,m=sin0sin#, andft=cos0. 
These substitutions give 



and 



A = c n sin 2 cos 2 <j>-{-- (c n ~ cu) sin 2 sin 2 #+ 

c 44 cos 2 0-t-2c i4 sin0 cos0 sin0, 
B = 7}(cu — c l2 ) sin 2 cos 2 + Cn sin 2 sin 2 

+ C44 cos 2 0— 2ci4 sin0 cos0 sine/), 
C=c 44 sin 2 + c 3 3 cos 2 0, 

F = c u sin 2 0(l — 2 sin 2 </>) + (C13+C44) sin0 cos0 sin<£, 
G=2cu sin 2 sin<£ cos$-f (C13+C44) sin0 cos0 cos</>, 

H = -(c u -\-c l2 ) sin 2 siii0 cos<£ + 2ci4 sin0 cos0 cos0. 



(22) 

(23) 
(24) 

(25) 
(26) 

(27) 
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To solve for the velocity of sound in a given direc- 
tion, the direction cosines must be substituted into 
the expressions for A ... H which must then be 
substituted into (15) and the resulting cubic equa- 
tion solved. For some directions this calculation is 
easy. For example, waves propagated along the 3 
axis will have /=0, m=0, n=l. 

Equation (15) reduces to 



c u — x 
c u — x 
c 33 — x 



which gives three roots 



pV — C33, C44, C44. 



(28) 



(29) 



These roots correspond to the propagation of a longi- 
tudinal wave with velocity V<Wp an( l two transverse 
waves with velocity ^c u /p. For other directions eq 
(15) does not reduce to such a simple form and the 
waves are quasi-longitudinal and quasi-transverse 
rather than truly longitudinal and transverse. 

A simplified method of treating directions for 
which eq (15) is difficult to solve has been used by 
each of the previous workers on the elastic constants 
of sapphire. The method has been most fully de- 
scribed by Mayer and Hiedemann [3]. They assert 
that for any direction the velocity of the quasi- 
longitudinal and quasi-transverse waves are given 
respectively by 

At>i=4a (30) 

and 

--cu (31) 



pvh 



where C33 and c' 4i are the value of c 33 and c u when 
the coordinate system is transformed into one with 
the X3 axis in the direction of propagation. Expres- 
sions for C33 and C44 are given by Mayer and Hiede- 
mann [3]. The present writers do not question the 
transformation equations for C33 and c' u but do 
assert that eqs (30) and (31) are not true in general. 
That is, pv\ or pv% will in general be a more compli- 
cated function of the elastic constants and direc- 
tion of propagation which may reduce to eqs (30) 
and (31) for certain directions of propagation. 

An example of a direction for which the simplified 
method gives a different answer than the general eq 
(15) is the Xi direction (parallel to a 2-fold axis) in 
sapphire. Using eqs (30) and (31) Mayer and 
Hiedemann obtain for pv 2 the three values 



and 



pvl=c n , 




(32) 


pv%=c u , 




(33) 


pV%=C^ = \{c n 


-Cu) 


(34) 



Actually substituting into eq (15) gives 



pvi=c lu 



pv%= 



B + O ± 1(^66— C 44 )X VH 



J 



(^66 C44) 



(35) 



(36) 



The results agree for the longitudinal wave but do 
not agree for the transverse wave unless c u is zero. 
In practice it might happen that 4ci 4 /(c 66 — c 44 ) 2 is 
small compared with one. If this were so, eqs (30) 
and (31) would give a good approximation for this 
particular direction but might not for a different 
direction. It seems that the results calculated from 

(30) and (31) would have to be checked against the 
results obtained from the fundamental eq (15), for 
each direction of propagation. This fact, of course, 
would mean that no time would be saved by the use 
of the simplified method. In the present work, only 
the exact eq (15) was used. 

It is perhaps worth noting why the simplified 
method should not be expected to be correct. If 
the XiX 2 x s coordinate system is transformed into the 
x[x' 2 x^ system, every term in eq (15) can he trans- 
formed to give a correct result. Equation (28), 
however, is a special case of (15) resulting from a 
special choice of coordinate system and direction of 
propagation. If the nonzero terms are transformed 
by a rotation of the coordinate system, eqs (30) and 

(31) will result, but these equations cannot be ex- 
pected to hold in general for directions of propagation 
other than the direction assumed in deriving eq (28). 

The error in the simplified method can perhaps be 
more easily seen from a physical argument. The 
transformed constants C33, C44, and c h r^ would be appro- 
priate for calculation of the velocities of waves 
traveling in the x' 6 direction and with particle motion 
along the x' s , x' 2 , and x[ axes, respectively. However, 
as Cady [16] has shown, these are not necessarily the 
directions of motion of the particles, and the devia- 
tion in some cases can be quite large. 

A method of solving eq (15) for the elastic con- 
stants which utilizes tensor transformations has been 
developed by Neighbours [28]. In his method all 
terms are taken into account, and in place of eqs 
(30) and (31) a set of equations is obtained each of 
which contains an infinite series which presumably 
converges rapidly. Equations (30) and (31) are 
equivalent to the zeroth approximation in this per- 
turbation method of Neighbours. He shows how 
rapidly the approximations converge in some cases 
but does not discuss the trigonal system. 

In his paper on calcite, Bhimasenachar [27] uses 
equations equivalent to (30) and (31) but explicitly 
states that small coupling terms have been neglected. 
In his subsequent paper on sapphire, Bhimasenachar 
[2] uses the simplified method without mentioning 
that it is approximate. Both Sundara Rao [1] and 
Mayer and Hiedemann [3] also use the simplified 
method but present it as if it were exact. 
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6. Specimens and Experimental Procedure 



specimens were of two 
ie dimensions, mass, 



The single crystal sapphire 

types: rods and blocks. T 

surface condition, and orientation of the 29 rod 
specimens are described in table 2. Two rectangu- 
lar block specimens were ground to our order by the 
Linde Company. The first block, SBl , was supposed 
to be 1 x}ix}i in. with tbe Long axis parallel to the x 3 
axis and one face perpendicular to the r, axis. 'Flic 
second block, SB2, was specified as lx%x% in. with 
the long axis in the .r 2 x^ plane making an angle of 135° 
with .7*3 and with one face perpendicular to Xi. The 
actual orientations are shown in figure 7. It will he 
seen that for SBl the 9 values are correct, but <j> is 6° 
instead of zero. For SB2 the (/> values are correct, 
but B is 47.5 instead of 45.0°. 



Table 2. Dimensions and orientation of specimens 


Specimen 


Mass a 


Nominal 
diametei 


Actual 
length 


Surface 

condition 


e 


4> 


SR 1 



2. 9870 
3. 1008 
3. 0737 
3. 0767 
2.1516 

3. 0647 
3. 0872 
3. 0902 
3.0636 

1.9297 

6. 9436 
6. 9374 
6. 9610 

6.9192 
6. 9099 

6.9151 
C). 9334 
6. 9219 
6. 9231 
4. 7551 

4. 7527 
3.1152 
3.1115 
3. 1024 
3. 1127 

3. 1154 
3.1155 
3. 1121 
3. 1173 


in. 

0.100 
.100 
. 100 
.100 
. 100 

. 100 

. 101) 

.100 
.100 
.100 

.150 
.150 
.150 
.150 
.150 

.150 
. 160 
.150 
.150 
.125 

.125 
.100 

. 100 
. 100 
. 100 

. 100 
.100 
.100 
.100 


cm 

14. 700 

15. 244 
15. 224 
15.232 
10. 664 

15.188 

15.218 
15.258 
15.240 
9. 591 

15.240 

i ;, -js2 

I5.2SO 
15.268 
15.288 

15.212 
15. 244 
15.268 

15. 280 
15. 260 

15. 228 
15.292 

15.2S2 
15.246 
15. 300 

15.296 
15. 300 
15. 296 
15.316 


f.p 

f.p 

f.p - 

f.p - 

f.p 

f.p 

f.p 

f.p 

f.p - 

f.p 

f.p 

C.g 

C.g 

C.g 

C.g 

C.g 

C.g 

C.g 

C.g 

C.g 

C.g 

f.p 

f.p 

f.p 

f.p 

f.p 

f.p 

f.p 

f.p 


13.6° 
10.7 
65.4 
29.4 
45. 6 

41.2 
67.8 
12.5 
89.0 
88.6 

59.5 
30 2 
11.9 
3.0 
18.8 

45.0 
40.7 
59.6 
64.0 
43.5 

48.8 
89.0 
89.5 
89.5 
89.0 

90.0 

ss.s 
89.5 
89.0 


+ 13.0° 


SR2 


+17.0 


SR3 


-2.8 


SK 1 


-24.0 


SR5 


+2.0 


SR6 

SR7 

SR8 


+2.0 

-8.0 

+20. 


SR9 


-0. 5 


SR 10 


+2.1 


SR 11 


-2.0 


SR 12 


+23.0 


SR13 


+27.0 


SR 14 


-30.0 


SR 15 ... 


+ 11.0 


SR16 -..- 


-29.0 


SR 17 


-28.3 


SR 18 


-28.0 


SR 19 


-1.0 


SR20 


-1.0 


SR21 


-0.5 


SR22 


+3.0 


SR23 


-3.5 


SR24 


-2.8 


SR25 


-3. 6 


SR26 


-3.7 


SR 27. . 


-2.5 


SR28_._. 


-3.0 


SR29.__. 


-3.0 







a Corrected for air buoyancy. 

Note: c.g. means centerless-ground. f.p. means flame polished after center- 
less grinding. The nominal diameters were not used in calculations. Instead, 
the average radius was computed from the actual length, mass, and density, 
P=3.9860±0.0010g/cm3. 



A Laue camera of convenl tonal design was used to 
determine the orientation; 35KVp X-radiation (cop- 
per target) was used with exposure times of about 
6 to 10 hr. The samples were supported with the 
specimen axis perpendicular to the X-ray beam and 
parallel to the edge of a small shield on the film 
holder. The shadow cast by this shield on the film 
established the reference line on the film. A speci- 
men to film distance of 3 cm was used, and the stereo- 
graphic projection was made according to the pro- 
cedure described by Barrett [29] using the Greninger 
chart given in his figure 22 on page 170, and a 21 -in. 
Wulff net. The resulting projection was then in- 




SBI 



SB2 



Figure 7. The orientation of the hvo rectangular block 
specimens. 



dexed with the letter symbols shown in figure 2a. 
This indexing was done 1 by a trial and error method 
of making a tentative assignment and comparing 
angles with the values for sapphire given by WincheD 
[8]. This process was continued until the proper 
fit was obtained. 

The second method of measuring 0, described in 
seel ion 2, was used. The actual construction obtained 
from a back reflection Laue pattern always resembled 
the projection shown in figure 8, that is, the opposite 
ends of the specimen axis were at opposite points of the 
stereographic equator. The value of d can be read 
directly from this projection, as shown, but the 
specimen axis must be projected on the X\X 2 plane to 
measure </>. Projection of +SA would require rota- 
tion of the whole pattern because PSA would lie on 
the hemisphere not shown in figure 8. This projec- 
tion can be done but requires more graphical con- 
struction with attendant possibility of error. It 
is more convenient to project — SA and choose X\ as 
the ± a direction nearest —PSA. The sign conven- 
tion for determining <f> is shown in figure 5. Orienta- 
tions were determined twice on 12 specimens with the 
resulting standard deviations A0 = O.()° and A<j>=\.Q°. 

The longitudinal resonance frequency of each rod 
shaped specimen was measured by hanging the speci- 
men vertically from a crystal pickup. A single, fine 
cotton thread was tied around the specimen some- 
what off center so that the specimen hung approxi- 
mately vertically. The other end of the thread was 
tied to the needle of a high output Eochelle salt 
pickup. Vibration was excited by placing a high 
frequency speaker just under the specimen and driv- 
ing through air. The electronic components used to 
( 1 ri \ e the speaker, to amplify t he pickup signal, and to 
measure the frequency have been described [30]. 
With this system it was possible to detect resonances 
in the range 30 to 50 kc. The response was weak 
but very sharp, and careful hunting was sometimes 
necessary. 



221 



-S A 




+ S A 



Figure 8. Stereogram illustrating the actual measurement of 
<f> and 0. 



The flexural resonance frequency was measured 
with the specimen suspended horizontally on two 
threads near the nodes. One thread was driven by a 
magnetic cutting head. The other thread was con- 
nected to a crystal pickup. Each specimen would 
vibrate in two flexural modes at right angles to each 
other. By turning the specimen in the string it was 
possible to excite either mode at will and to determine 
the direction of vibration for each. The two 
resonance frequencies usually differed by a few 
cycles, presumably because of small variations in the 
diameter of the specimen. The average of the two 
resonance frequencies was used in eq (13). An 
average value of r 2 , for each rod calculated from the 
density, the mass, and the length, was used. This 
approach was used because the density could be 
determined accurately by hydrostatic weighing, but 
the average radius was more difficult to measure 
accurately. 

The torsional resonance frequency was determined 
with the same speaker and pickup used in the longi- 
tudinal resonance frequency measurements. The 
tweeter was modified by covering its output opening 
with a sheet of paper. The thread used to drive 
the specimen was fastened to the center of this sheet. 
The thread was wrapped around the specimen several 
times to provide an off-center driving force and the 
pickup thread was wrapped around the specimen in 
the opposite direction. Careful hunting was again 
necessary to find the weak but sharp resonance. 
It was necessary to distinguish between the funda- 
mental of torsion and the overtones of flexure. This 
distinction was accomplished by moving the driving 
and pickup strings until the nodes were located. 
The counter used for frequency measurements had an 
accuracy of ±0.1 cps. 

The velocity of sound measurements on the two 
block specimens were made by sending an ultrasonic 
pulse into the specimen and measuring the transit 
time. A quartz crystal 0.5 in. in diameter and cut to 
resonate at 10 Mc was used as transmitter and pick- 



up. An X-cut crystal was used to generate longi- 
tudinal waves, and an AC-cut crystal was used for 
transverse waves. These crystal cuts are described 
by Buchanan [31]. The crystals were mounted on 
the specimen with salol (phenyl salicylate). The 
specimen was coated with a very thin layer of evapo- 
rated aluminum before the crystals were mounted in 
order to provide a ground. 

The crystals were excited and the echoes were am- 
plified with a commercial ultrasonic generator and 
receiver. This unit had a pulsed oscillator of variable 
frequency with a fixed pulse repetition rate of 100/ 
sec. The amplified output of the oscillator was con- 
nected to an electrode on a quartz crystal. The same 
electrode served as an input connector for the re- 
ceiver. The output of the receiver was connected to 
one input of an oscilloscope with a dual trace pre- 
amplifier. A crystal controlled 1-Mc oscillator was 
connected to the other input of the dual trace pre- 
amplifier. The oscilloscope thus displayed the stand- 
ard 1-Mc signal on one trace followed by the echo 
pattern on the next trace. These two traces were 
simultaneously visible, and the 1-Mc pattern was 
used as a time scale for the transit time measurement. 
The general technique of pulse velocity measurement 
is reviewed by Huntington [32]. 

The sweep expansion feature of the oscilloscope 
could be read to 0.01 jji sec. The actual measurements 
were much less accurate, however, because the pulses 
were rounded and did not have a sharply defined lead- 
ing edge. Consequently the maximum of the pulse 
envelope was used as a reference point. Repeated 
measurements suggested that the resulting transit 
time might have an error of several percent. For 
this reason no transit time correction, such as that 
discussed by Overton and Gaffney [33], was made. 
The measured velocities given in this report should 
be examined with this in mind. They are provided 
as a rough check on the accurate values determined 
by the resonance methods to show that the values 
of the elastic constants, within the accuracy of 
pulse velocity measurements, show no frequency 
dependence. 

7. Results 

The resonance frequencies are summarized in 
table 3. These values were used to calculate E f and 
G v using the equations of section 4. The method de- 
scribed in section 3 was then used to calculate the 
elastic compliances, s tj . Tables 3 through 7 show 
various stages in these calculations and indicate the 
precision. The reciprocals of the E f values used in 
solving eq (8) are given in table 4, and the resulting 
values of the compliance combinations are given in 
the first half of table 6. The difference column in 
table 4 shows how well the observations on all 29 
rods were fitted by T the four parameters of eq (8). 
The shear modulus calculations are summarized in 
table 5, and the resulting compliance combinations 
are given in the second half of table 6. The difference 
column of table 5 shows how accurately eq (9) fitted 
the data. 
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Table 3. 


Resonance frequencies 


of sapphire 


rods 


Specimen 


Flexural :i 


Longitudinal 


Torsional 




/i 


ft 




SRI... 


1100. l 
L032.fi 
971.fi 

990. 
1908. 9 

945.1 
998. 3 
1021. 9 

993.8 
2518.0 

1431.9 

1378. 7 
1536.0 
1557. 1 
1487. 6 

1483.0 

1472. 5 
1508.0 
1433.7 
1166.7 

1172.7 
1015.4 
1015.9 

1022. 6 
1015.0 

1013. 5 
1007. 7 
1013. 1 
1010.6 


1100.7 
1033. 7 
972; 6 
980.9 

1907.3 

943.6 

997. 5 

1023.7 

1010. 3 

2fi39. 3 

1432. 7 
1380. 1 
1537.0 


35673 
34739 
32747 
33397 

45060 

31746 
33538 
34464 
33868 

53817 

32200 
31145 
34607 
35175 
33703 

33336 
33160 
34068 
32420 
31778 

31819 
33727 
33786 
33871 
33768 

33735 
33757 
33715 
33790 


20850 


SR2... 


19940 


SR3... 


21033 


SR4... 


20547 


SR5... 




SR6 


21600 


SR7... 


20616 


SK s 


20037 


SK 9 


20505 


SR 10— 




SR 11 


21210 


SR 12 


21763 


SK 13 


19943 


SR14 


19704 


SR 15... 


1488.4 

1484.1 
1473. 5 
1507.7 
1434. 1 
1169.2 

1174.9 
989.0 
992. 2 
999. l 

989. 9 

98ft 2 
99a 

986. 3 
995. 6 


20346 


SR16 


20596 


SR 17 


20587 


SK IS 


20152 


SR 19 


20992 


SR20 


21331 


SB 21... 


21436 


SR22... 


20251 


SR23. .. 


20247 


SR24... 


20332 


SR25... 


20236 


SR 26... 


20230 


SR27... 


20219 


SR28._- 


20233 


8B29... 


20240 







:i The two flexural frequencies correspond to vibrations at right angles. The 
average value was used for calculations. All values are in cycles per second. 



Table 4. Reciprocal Young's modulus of sapphire 

All values are in units of 10~ 13 cm^dyne. The standard deviation for 1/E, 
computed using the tabulated differences is ±0.0052X10~ 13 cm 2 /dyne. 



Specimen 


1/Ef 
Flexural 


1/Ef 

Longitudinal 


1/E/ 

Calculated •■ 


Difference * 


SRI: 

SR2 


2. 2764 
2. 2392 
2. 5242 
2. 4262 
2. 7124 

2. 6968 
2. 4088 
2. 2684 
2. 3530 
2. 3543 

2. 6020 
2. 7684 
2. 2505 
2. 1795 
2. 3679 

2. 4414 
2. 4570 
2. 3214 
2. 5564 
2. 6689 

2. 6709 
2. 3537 
2. 3488 
2. 3424 
2. 3440 

2. 3535 
2. 3482 
2. 3436 
2. 3486 


2. 2805 
2. 2366 
2. 5233 
2. 4237 
2. 7167 

2. 6976 
2. 4079 
2. 2681 
2. 3540 
2. 3540 

2.6048 
2. 7685 
2. 2461 
2. 1744 
2. 3624 

2. 4390 
2. 4546 
2. 3180 
2. 5556 
2. 6674 

2. 6717 
2. 3579 
2. 3527 
2. 3521 
2. 3497 

2. 3554 
2. 3512 
2. 3477 
2. 3521 


2. 2759 
2. 2369 
2. 5193 
2. 4277 
2. 7119 

2. 6876 
2. 4147 
2. 2621 
2. 3529 
2. 3564 

2. 5967 
2. 7708 
2. 2545 
2. 1744 
2. 3632 

2. 4404 
2. 4567 
2. 3125 
2. .5645 
2. 6670 

2. 6826 
2. 3560 
2. 3515 
2. 3518 
2.3502 

2. 3529 
2. 3509 
2.3517 
2.3507 


0. 0046 
-.0003 

. 0040 
-. 0040 

.0048 

.0100 

-.0068 

. 0060 

.0011 

-.0024 

.0081 
-.0023 
-.0084 

. 0000 
-.0008 

-.0014 
-.0021 

.0055 
-.0089 

. 0004 

-.0109 
.0019 
. 0012 
.0003 

-.0005 

.0025 

.0003 

-.0040 

.0014 


SR3 

SR4 

SR5 

SR6 

SR7 

SR8 

SR9 

SR 10. 


SR11 


SR 12 


SR 13 


SR 14 


SR 15. . 


SR 16 


SR 17 


SR 18 


SR 19 


SR20 


SR21 


SR22 


SR23 


SR24 


SR25 


SR26 


SR27 


SR28-. 


SR29 





° Calculated from eq (7) using the data of table 6. 

b Difference= longitudinal value minus calculated value. 
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Table 5. Reciprocal shear modulus of sapphire 

All values are in units of 10~ 13 em 2 /dyne. The. standard deviation for 1/0/ 
calculated from the dilTerencc column is ±.()178X10- 13 cm 2 /dyne. 



Specimen 


(««)*+ (*m) 2 


1/0/ 
Torsional 


1/Qf 
Calculated 


Difference 


2 * S 33 


SR 1 


0. 0421 
. 0291 
,0716 
.0124 


6.7186 
''..SI 75 
6.1888 
6.4166 


6. 7296 
6. 8066 
6. 2242 
6. 4294 


—0 0110 


SR 2_._ 


0109 


SR 3... 


— 0354 


SK 1 


— 0138 


SR 5» 




SR 6 


.0521 
. 0564 
. 0390 
. 0002 


5. 8792 

6. 4282 
6. 7491 
6. 5208 


5.9099 
(i. 4295 
0. 7569 
6. 5396 


— 0307 


SR 7 


— 0013 


SR 8. . 


— 0078 


SR 9 _. 


— 0188 


SR 10* 




SR 11 


.0734 
. 0064 
.0371 
. 0020 
.0643 

.0049 
. 0010 
. 0189 
.0748 
. 0506 

. 0564 
. 0024 
. 0014 
. 0008 
.0010 

.0019 

. 0004 
. 0009 
.0007 


C). 0765 

5. 7365 

6. 8003 
6. 9320 
6. 5469 

6. 3939 
6. 3692 
6.6439 

6. 1704 

5. 969S 

5. 9428 

6. 5421 
6. 5528 
6. 5282 
6. 5439 

6. 5520 
C). 5544 
6. 5445 
6. 5319 


6.0763 
5. 7468 
6. 7718 
6. 9302 
6. 5571 

6. 3990 

(). 3690 
6. 6407 
6. 1355 
5. 9498 

5. 9153 
6. 5334 
6. 5425 
6.5419 
6. 5451 

6. 5396 
6. 5435 
6. 5421 
6. 5440 


0002 


SR 12 ... 


— 0103 


SR 13 


0285 


SR 14 


0018 


SR 15 


— 0102 


SR 16 

SR 17.. . 


-.0051 
0002 


SR 18 

SR 19 


.0032 
0349 


SR 20 


0200 


SR 21 

SR 22 


.0275 

. 0087 


SB 23 


0103 


SR 24.. .. 


— 0137 


SR 25..-. 


— 0012 


SR 26 


0124 


SR 27- 


0109 


SR 28. 


0024 


SR 29 


— 0121 







<■ Too short for use. 



Table 6. Best estimates of some 


compliances 


Compliance 


Value 


Standard deviation 


Calculated from Young's modulus values <» 


*33 

*44+2*l8 

8U 


2.3529X10- 13 cm-'/dvne 
2. 1694 
6. 2183 
0. 4901 


0.0016X10- |3 cm-7dvne 
. 0024 
. 0088 
.0057 


Calculated from shear modulus values * 


8n— «i2 

944 

311+533- 
814 


(2St3-H44) 


3.0696X10-'3cm2/dyne 
6. 9400 
-1.6539 

0.4868 


0.0067X10-' 3 cm-'/dyne 

.(M)Nl 

.0197 
.0101 



u Calculated by fitting eq (7) to Young's modulus values calculated from 
longitudinal resonance frequencies. 

b Calculated by fitting eq (8) to shear modulus values. 



The final result for the elastic compliances and 
constants is given in table 7. The six compliances 
were calculated from the eight values in table 6 by 
weighting each value according to its standard de- 
viation. The constants were calculated from the 
compliances using the standard matrix inversion 
equations given, for example, by Nye [7]. 

The Cfj values of table 7 were used in eq (15) to 
calculate the velocity of sound for comparison with 
the measured velocities for all directions of propa- 
gation normal to the faces of SBl and SB2. In 
some cases this necessitated the solution of a cubic 
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equation by successive approximations. The re- 
sults are given in table 8 along with the measured 
values. In this table the faces have been described 
by the labels shown in figure 9, e.g., the 45YZ face. 
The actual orientation was used in the calculations, 
however; that is, the 45 YZ face normal actually had 
61=47.50° instead of 45°. 

The Sij values were used to calculate E f and G f as a 
function of orientation. The results are shown in 
figures 9 and 10 and in contour form in figures 11 
and 12. 

8. Discussion 

The values of the elastic constants and compliances 
for sapphire which have been reported by other 
investigators are given in table 7. Mayer and 
Hiedemann have suggested that the method of 

Table 7. Elastic compliances and constants of sapphire 
All sn values in units of 10- 13 cm 2 /dyne. All c,/ values in units of 10 12 dyne/cm2. 





Sundara 
Rao 


Bhimase- 
nachar 


Mayer & 
Hiedemann 


Present work 


»1I___ __ 


2.84 
2.21 
5.47 

-0.95 
-.47 

-1.52 

4.66 
5.06 
2.35 
1.27 

1.17 
0.94 


2.32 
1.93 

5. 77 
-1.05 
-0.38 
-1.71 

4. 65 

5.63 
2.33 
1.24 
1.17 
1.01 


2.18 
2.02 
5.04 
-0. 50 
-. 16 
-.49 

4.96 
5.02 
2.06 
1.09 
0.48 
.38 


2. 353±0. 002 




2. 170± . 002 


«44 


6.940db .008 


Si2 


-0.716± .007 




-.364± .006 




.489± .005 


C\\ . ... 


4. 968±0. 018 




4.981± .014 


C44 -- --- - 


1.474± .002 


Ci2 - - - - 


1.636± .018 


Cm . _ _ 


1. 109± .022 


C\4 -- 


-0.235+ .003 







Table 8. Velocity of sound in sapphire 

All velocities in 10 3 cm/sec. 



Direction ! 



Mode b 



Calculated 
velocity 



Measured 
velocity 



Percent 
difference 





Specimen SB1 






X 


L 
Tpol Y 
Tpol Z 

L 
T pol A^ 
Tpol Z 

L 
T 


1117 
667 

585 

1118 
651 
600 

1118 
608 


1092 
669 
579 

1116 

645 
609 

1100 
594 


-2.2 


Y ... 


+0.3 
-1.0 

—0.2 


Z 


-0.9 
+1.5 

— 1.7 




-2.3 



Specimen SB2 



45° YZ.. 



135° YZ. 



L 

Tpol Y 
Tpol Z 

L 

Tpol X 

T other 

pol 

L 

T pol A' 

T other 

pol 



1116 
676 
575 

1C95 
581 
691 



1041 
671 



662 
564 

1146 
609 
710 



1040 
659 
687 



-2.5 
-2.1 
-1.9 

+4.7 
+4.0 
+2.7 



-0.1 
-1.8 
-0.4 



a The present directions are the actual face normals shown in figure 8. 

b ".L" means longitudinal, " T pol Y" means transverse with the direction of 
vibration along Y, and " T other pol" means transverse with the direction of 
vibration perpendicular to the other two modes. 
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Young's modulus, Ef, as a function of orientation 
for single crystal sapphire. 
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Figure 10. 



The shear modulus, Gf, as a function of orientation 
for single crystal sapphire. 
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exciting vibrations used by the Indian workers may 
give spurious resonances. For this reason we com- 
pare primarily with the set of values given by Mayer 
and Wiedemann. 

The met hod used by Mayer and Hiedemann gave 
the Cij values directly; an error in a single one of the 
Cij might lead to errors in several of the *,,■ values. 
Accordingly we compare c values rather than s 
values. c u and c 83 are in good agreement, but the 
other values are not. The difference in sign of r,., 
is striking. As was pointed out earlier, the sign of 
c u (and s H ) depends on which end of the 2-fold axis 
is taken as -f,/:,. The conventions given by Xye |7] 
and by the Standards on Piezoelectric Crystals [11] 
leave this unspecified. This point is not mentioned 
by any of the previous workers on sapphire, and we 
have not been able to find any information in their 
papers on which choice was made. If Mayer and 
Hiedemann used the opposite end of a 2-fold axis as 
-\-x ly the difference in sign of c 14 would be accounted 
for and the values of Cu would be in fair agreement. 
The most troublesome disagreement is thai between 
the value of r,, obtained in the present work, 1.47X 



10 12 dyne/cm 2 , and the value of 2.06X10 12 dyne/cm 2 
obtained by Maver and Hiedemann. The 1.47 value 
cannot be very much in error because the value of 
the velocity of the transverse wave in the x 3 direction 
depends only on this constant, and the measured 
value is —2.3 percent from the calculated value as 
shown in table 8. The fact that eqs (30) and (31), 
used by Mayer and Hiedemann, are not strictly cor- 
rect was pointed out in section 5. However, it does 
not seem that this fact has any connection with the 
disagreement in the values of c 44 because these 
equations are not needed Fora wave propagated in the 
£ 3 direct ion. ' 

A comparison can be made to support the validity 
of the present data. Lang [34] has reported values 
of elastic moduli determmed dynamically on hot 
pressed polycrvstalline alpha alumina of density 
3.983 g/cm 3 . This density is so close to the 4 value of 
3.986 g/cm 8 obtained on single crystals that the 



itin 



i dated March 28, I960, R. H. Truell (Brown 
wits from his joint Investigation with P. H. Levy 
^iwjuiqvcu .>,...«Miai i «*■ »wi at ory) using a pulse velocity technique, C33= 5.013 
Xio 1 - dynes/cm* and C44=1.478X10 ,a eynes/cm 2 . These values are seen to be in 
good agreement with the results presented here. 



4 In a private 
Universit3 1 givi , .. 
(Brookhaven Nat tonal Lab 




Figure 11. Young's modulus, K f , intervals of 100 kilobars. 
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values of the elastic moduli of his specimens must 
be only slightly less than the values for zero porosity. 
He gives £'=4012 kilobars and 6r=1591 kilobars as 
his best values. An exact theoretical relationship 
between single crystal elastic constants and poly- 
crystalline moduli has not been developed, but two 
good approximations are discussed by Huntington 
[32]. One method, due to Voigt, should give values 
slightly too high, and the other method, due to Reuss, 
should give values slightly too low. Using the elastic 
constant and compliance data of table 7 and the 
equations given in Huntington's article, one obtains 
the results shown in table 9. The values of E and G 
calculated from the elastic constants obtained in this 
investigation are in rather good agreement with 
Lang's values. 

All elastic constants, compliances, and moduli so 
far discussed have been for adiabatic conditions; i.e., 
measured under conditions when appreciable heat 



flow could not occur during one cycle of vibration. 
Isothermal values are measured in static tests. 
Hearmon [14] gives an equation from which the differ- 
ence between isothermal and adiabatic values can 
be calculated. This difference ranges from 0.2 
percent for s n to 1 percent for s 13 . 



Table 9. 



Comparison of calculated and measured elastic 
moduli jot poly crystalline alumina 



All values in kilobars. 



Measured (Lang) 

Present 
data 

Mayer & 
Hiedemann 



Reuss theory. 
Voigt theory. 

Reuss theory. 
Voigt theory 



Young's 
modulus 



3973 
4084 



4554 

4682 



Shear 
modulus 



1606 
1660 



1996 
2070 



1440 




Figure 12. Shear modulus, Gf, intervals of 20 kilobars. 
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The compressibilities can be calculated from the 
compliances using equations given by Nye [7]. The 
adiabatic compressibilities calculated in this way are 
given in table 10. These should differ from the 
isothermal compressibilities by about 1 percent. 
Values of the isothermal compressibilities measured 
in static tests by Bridgman [35] and by Madelung 
and Fuehs [36] are given in table 10. Bridgman's 
linear compressibility parallel to the x z axis is based 
on the initial slope of his curve and is inconsistent 
with all the other data. It is listed in the first row 
of measured values in table JO even though it does 
not appear to be consistent with his reported value 
of volume compressibility, [f the slope of the sub- 
sequent section of Bridgman 's curve 4 is used, the 
values labeled "Bridgman, interpreted" are obtained. 
The present work is in somewhat better agreement 
with the experimental values of volume compress- 
ibility than the results of Mayer and Hiedemann. 
Visual inspection of Bridgman' s curves suggests that 
the difference between the value of 3.6X10" 13 cm 2 / 
dyne, calculated from his data, and 3.988X10~ 13 
cm 2 /dyne may be largely due to experimental error 
in his data. 



Table 10. Comparison of compressibilities calculated from 
elastic compliances with values measured directly 



All values in LO-is cm 2 dyne. 





Linear compress- 
ibility 


Volume 




Parallel 
to.r 3 


Perpen- 
dicular 

tO.T3 


compress- 
ibility 




Calculated from compliances 


Present data, 25° C 


1.442 
1.70 


1.273 

1.52 


3.988 


Mayer and Hiedemann, 27° C._ 


4.74 




Measured 


Bridgman, 30° C, as reported... 
Bridgman, 30° C, interpreted... 
Madelung and Fuchs, 0° C 


0. 40 
1.35 


1.13 
1.13 


3.23 
3.6 

3.8 











Finally, the difference columns in table 4 and table 
5 show how well the classical, 6-constant theory 
fits the accurate results obtained in the kilocycle 
per second range. The standard deviation calculated 
from either set of differences is 0.3 percent of the 
average value of \\E or 1/6?, respectively. The 
worst disagreement for either \\E or 1/6? is only 
0.6 percent. It is evident that the classical theory 
of elasticity gives a very accurate representation of 
the elastic behavior in the kilocycle per second range. 



9. Summary 

(1) The complete set of elastic constants and 
compliances for sapphire calculated from resonance 
measurements in the kilocycle per second frequency 
range is given. 



(2) Less accurate pulse velocity measurements in 
megacycle per second range indicate little, if any, 
frequency dependence of the elastic constants. 

(3) The values of r M and r :vi agree well with the 
results of the most accurate previous measurements, 
hut the remaining constants and all of the compli- 
ances are in serious disagreement. 

(4) Values of Young's modulus and the shear 
modulus for polvcrvstalline alumina calculated from 
the compliances are in good agreement with experi- 
mental values. 

(5) Curves for estimating Young's modulus and 
the shear modulus for single crystal sapphire of any 
orientation are given. 
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